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Abstract. The nature of this paper is twofold: On one hand, we will give a short introduc- 
tion and overview of the theory of model sets in connection with nonperiodic substitution 
tilings and generalized Rauzy fractals. On the other hand, we will construct certain Rauzy 
fractals and a certain substitution tiling with interesting properties, and we will use a new 
approach to prove rigorously that the latter one arises from a model set. The proof will 
use a duality principle which will be described in detail for this example. This duality is 
mentioned as early as 1997 in |Gel| in the context of iterated function systems, but it seems 
to appear nowhere else in connection with model sets. 



1. Introduction 

One of the essential observations in the theory of nonperiodic, but highly ordered structures 
is the fact that, in many cases, they can be generated by a certain projection from a higher 
dimensional (periodic) point lattice. This is true for the well-known Penrose tilings (cf. fGSH') 
as well as for a lot of other substitution tilings, where most known examples are living in , 
or E^. The appropriate framework for this arises from the theory of model sets ( |MEYj . 
IMooj ). Though the knowledge in this field has grown considerably in the last decade, in 
many cases it is still hard to prove rigorously that a given nonperiodic structure indeed is a 
model set. 

This paper is organized as follows. Section [l.ll briefiv collects some well-known facts about 
substitution tilings. These facts will be used later. Section 11.21 contains some basics about 
model sets, in particular in connection with substitution tilings. In Section 1.3 we give the 
definition of Rauzy fractals. This can be done completely in the framework of model sets; it 
is just a question of the terminology. In Section we introduce a family of one-dimensional 
substitutions which are generalizations of the well-known Fibonacci sequences. Four of these 
possess a substitution factor which is a PV-number, so these are candidates for being model 
sets. Section 12.21 is dedicated to one of these substitutions. Therein, the corresponding 
Rauzy fractal is obtained. Section [2.31 shows the construction of the dual substitution tiling, 
which is a two-dimensional tiling with some interesting properties. Using the duality, we 
will show that this dual tiling is a model set. Hence, this is one of the few known cases of 
a two-dimensional model set, generated by a substitution, where the substitution factor is 
not an algebraic number of degree one or degree two (and the proof being rigorous). Section 
01 contains some additional remarks, and further examples of dual substitutions, where the 
duality principle can be used in order to prove that these are model sets, too. 
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Let US fix some notation. Throughout the paper, cl(M) denotes the closure of a set M, int(M) 
the interior of M, dM the boundary of M, and #M the cardinality of M. denotes the 
d-dimensional Euclidean space, i.e., equipped with the Euclidean metric || • ||. B*^ denotes 
the unit ball {x \ \\x\\ < 1}. If not stated otherwise, /x denotes the d-dimensional Lebesgue 
measure, where d will be clear from the context. A tile is a nonempty compact set T C E'^ 
with the property cl(int(T)) = T. A tiling of E"' is a collection of tiles {r„}^>Q which covers 
E*^ and contains no overlapping tiles, i.e., int(T^) n int(T„) = for k ^ n. In other words: 
A tiling is both a covering and a packing of E'^. A tiling T is called nonperiodic, if the only 
solution of T + £C = T is £C = 0. 

1.1. Substitution Tilings. A convenient method to generate nonperiodic tilings is by a 
substitution: Choose a set of prototiles, that is, a set of tiles J- := {Ti,T2, . . . ,Tfn}. Choose 
a substitution factor A G M, A > 1 and a rule how to dissect ATj (for every 1 < i < m) into 
tiles, such that any of these tiles is congruent to some tile Tj G T . It will be sufficient for this 
paper to think of a substitution like in Fig. ^ For completeness, we give the proper definition 
here. 

Definition 1.1. Let T := {Ti,T2, . . . ,Tm} be a set of tiles, the prototiles, and X > 1 a real 
number, the substitution factor. Let XTi = (pi(Ti-^)L)ip2{Ti,^)L) . . .L)ipn{i)(Ti^^^-^), such that every 

tpj is an isometry of K'^ and the involved tiles do not overlap. Furthermore, let a{{Ti}) = 
{931 (Tjj), Lp2{Ti^), . . . , (Tj^j.j )}. Let S be the set of all configurations of tiles congruent to 
tiles in T . That is, S is the set of all sets of the form {Tj^ + Xi \ Tj- G T ,i G /, iCj G E"^}, 
where I is some finite or infinite index set. By the requirement cr({Tj + a;}) = a{{Ti\) + \x 
(x £ E'^J, a extends in a unique way to a well-defined map from S to S. Then a : S —>■ S is 
called a substitution. Every tiling of K'^, where any finite part of it is congruent to a subset 
of some a^{Ti) is called a substitution tiling (with substitution a). 

More general, this definition makes sense if we replace A by some expanding linear map. That 
is, replace A in the definition above by a matrix Q such that all eigenvalues of Q are larger 
than 1 in modulus. In this full generality, many of the following results may break down. 
Anyway, in Section [2. 21 we will use the generalization in one certain case, where it will cause 
no problem. 

To a substitution a we assign a substitution matrix: := {aij)i<ij<m, where Oij is the 
number of tiles Tj + x (tiles of type Tj G J^) in cr({Tj}). In this paper we will consider 
substitutions with primitive substitution matrices only. Recall: A nonnegative matrix A is 
called primitive, if some power A^ is a strictly positive matrix. In this situation we can apply 
the following well-known theorem. 

Theorem 1.2 (Perron-Frobenius) . Let A G K"^^™ 5e a primitive nonnegative matrix. Then 
A has an eigenvalue < A G M, which is simple and larger in modulus than every other eigen- 
value of A. This eigenvalue is called Perron-Frobenius-eigenvalue or shortly PF-eigenvalue. 
The corresponding eigenvector can be chosen such that it is strictly positive. Such an eigenvec- 
tor is called PF-eigenvector of A. No other eigenvalue of A has a strictly positive eigenvector. 

Using this theorem, it is a simple exercise to show the following facts (cf. |PYFj . Thm. 1.2.7): 
We consider the relative frequency of the tiles of type Tj in a tiling T, i.e., the ratio of the 
number of tiles of type Tj in T and the number of all tiles in T. Formally, this is defined by 
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. Since we consider primitive substitution tilings only, it is true 



oo 



#{T+x er I xdrB'i , Tar} 



that this limit exists. Moreover, this limit exists uniformly^ i.e., for any translate T + a; we 
will obtain the same limit. 

Proposition 1.3. Let a be a primitive substitution in E"^ with substitution factor A and pro- 
totiles Ti, . . . , Tm- Then the PF-eigenvalue of a Aa- is X'^. In particular, a substitution factor 
is always an algebraic number. 

The normalized (right) PF-eigenvector v = [vi, . . . of A(j contains the relative frequen- 

cies of the tiles of different types in the tiling in the following sense: The entry Vi is the 
relative frequency ofTi in T . 

The left PF-eigenvector (resp. the PF-eigenvector of ) contains the d-dimensional volumes 
of the different prototiles, up to scaling. 

Some people don't like the term 'left eigenvector'. Those may replace 'left eigenvector of A^ 
by 'eigenvector of A'^ , written as a row vector' wherever it occurs. In the next section, we 
will associate with a tiling a related point set: By replacing every prototile Tj with a point in 
its interior, a tiling T gives rise to a Delone set Vr- (Note that this defines Vr not uniquely. 
It depends on the distinct choice of points in the prototiles. But the following is true for 
any Delone set Vr constructed out of a primitive substitution tiling T in the described way.) 
A Delone set is a point set V <Z which is uniformly discrete (there exists r > such 
that #(1^ n rB'^ + a;) < 1 for all x G E"^) and relatively dense (there exists > such that 
#(y n m'^ + x)>l for ah a; G E'^). Since the frequencies of the prototiles tiles (cf. remark 
after Theorem II. 2|) are well-defined for all tilings T considered in this paper, it makes sense 

to define the density of the associated Delone set Vr'- dens(Vr) = limr^oo ^^^^fSy"- 1^ is not 
hard to show the following identity. 

Proposition 1.4. Let T be a substitution tiling with substitution a, let v be the right nor- 
malized PF-eigenvector of A^ (containing the relative frequencies of the prototiles), let tv^ 
be the left PF-eigenvector of A^, such that w'^ contains the d-dimensional volumes of the 
prototiles, and Vt a Delone set constructed out of T by replacing every prototile with a point 
in its interior. Then 



Example: The golden triangle substitution. The substitution rule, together with substitution 
factor, substitution matrix etc. is shown below. 



Figure 1. 'Golden triangles'. A substitution rule in E^ with two prototiles 
Ti , T2 and substitution factor A = where r = ^ is the golden mean. 

This substitution uses two prototiles Ti , T2 . The PF-eigenvalue of the substitution matrix Aa- 
is r = A^. The normalized PF-eigenvector is v := (r~^,r~^)"^. Thus, the relative frequency 
of the tiles of type Ti (resp. T2) in the tiling is (resp. r^^). A^ is symmetric, thus 
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v'^ is a left PF-eigenvector. The areas (2-dim volumes) of the two prototiles are /^(Ti) = 
-^/r^, /i(T2) = -y/r. If we replace every triangle in a corresponding substitution tiling by a point 

in its interior, then by Prop. IlH the density of the point set is (^{\/t^, \/t) ■ {t~^,t~'^^'^ 

T-3/2 



(7^ 



0.56886448 . . . 



1.2. Model Sets. Model sets are a special kind of Delone sets. They show a high degree 
of local and global order. For example, a model set V and its difference set V — V := 
{x — y\x,y £ V} do not differ too much in the following sense (cf. [Mey]): There is a finite 
set F such that 

V -V + F. 

Lattices do have this property with F = {0}, and indeed lattices are special cases of model 
sets. Unfortunately, there is no reference known to the author which gives a comprehensive 
overview about model sets. A good reference is |Mooj . A comprehensive overview for the 
case of one-dimensional tilings is contained in |PYFj . mainly chapter 7 and 8. Therein, many 
connections to combinatorics, ergodic theory and number theory are compiled (but the term 
'model set' cannot be found in this book, the term 'geometric representation' is used instead). 

A model set is defined by a collection of spaces and maps (a so-called cut-and-project scheme) 
as follows. In let A be a lattice of full rank in E*^^*^, 7ri,7r2 projections such that ttiIa 
is injective, and 7r2(A) is dense in E*^. Let be a nonempty compact set — the so-called 
window set — with the properties cl(int(VF)) = W and fi{dW) = 0, where fi is the Lebesgue 
measure in E*^. 

(1) u u u 

V A w 

Then the set V = {7ri(a;) | x G A, 7r2(a;) G W} is called a model set^. Obviously, y is a set of 
points in E*^. A tiling can be obtained by applying certain rules to V, e.g., considering the 
Voronoi cells of V. In the case d = 1 it is quite simple to get a tiling: If V = {xi \i £ Z}, 
such that Xi < ajj+i, then T = {[xi,Xi^i] | i G Z} is a tiling of E^. 

Since V C vri(A) and vtiIa is injective, every x £ V corresponds to exactly one point tt^^{x) G 
A, and thus exactly to one point 7r2(7r]^"'^(a;)) G W. Therefore, we are able to define the star 
map: 

(2) * : ^i(A) ^ E^ X* = 7T2iTT^Hx)) 

This allows for a convenient notation. E.g., we have 

(3) d{V*) = W, 
where V* = {x* \ x eV}. 

Remark: In general, E*^ and E*^ in can be replaced by some locally compact abelian groups 
G and H. Then ^ denotes the Haar measure of H. This is necessary to show that, for 
example, the chair tiling is a model set (cf p3 MScl ) . If the substitution matrix A^j of a given 
substitution tiling is unimodular (i.e., det(Ao- = if), then it suffices to consider Euclidean 
spaces only. If not, then one needs to consider more general groups (cf. [S])- In this paper, all 

^To be precise, this defines a regular model set. Without the property iJ.{dW) = 0, is still a model set, 
but not a regular one. 
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occurring substitution matrices are unimodular, therefore we will consider Euclidean spaces 
only. 

Given a substitution tiling, one may ask if this tiling can be obtained as a model set V. By 
the following theorem, this is only possible if the substitution factor is a PV number or a 
Salem number. 

Definition 1.5. Let A G M 6e an algebraic integer. A is called a Pisot Vijayaraghavan number 
or shortly PV number, if \X\ > 1, and for all its algebraic conjugates Aj holds: \Xi\ < 1. 
If for all conjugates Aj of X holds |Aj| < 1, with equality for at least one i, then X is called a 
Salem number. 

Theorem 1.6 (Meyer). If V is a model set, 1 < A G R and if XV C V , then X is a PV 
number or a Salem number. 

It is known that every substitution a gives rise to at least one tiling T with the property 
cj(T) = T. So, if F is a point set derived from T (where the substitution factor is A), say, 
V being the vertex set of this tiling, then by the definition of a substitution follows XV C V . 
Now, if V is also a model set, by the theorem above the substitution factor must be a PV 
number or a Salem number. 

The remaining things to check in order to show that F is a model set are the following (cf. (fUl): 

• Determine the lattice A, 

• determine the window set W , 

• show ii{W) > 0, 

• show n{dW) = 0, 

• show dens(y) = dens({7ri(ic) | x G A,tt2{x) G W}). 

There are standard constructions for obtaining A and W. But the latter three items are hard 
to decide rigorously in general, or even in special cases. For an example how hard it may be 
to give a proper proof, cf. |BSj . In the next section, we will use a duality principle to prove 
that the three latter items are fulfilled for the considered tiling, thus proving that it arises 
from a model set (i.e., the associated point set is a model set). But first, we will discuss the 
above list in more detail. 

Starting with a substitution tiling, or better a substitution point set V in (cf. the remark 
after (0)), the standard construction for the lattice is as follows: Let A be the substitution 
factor (where A from now on is always assumed to be a PV number) and X2,...,Xra its 
algebraic conjugates. Then 

^ = ((1) 1) ■ • • 1 1)^) (X, A2, • • • , Am)^, (A^, A2^, . . . , Am^)^, • • • , {X"^'^, X2'^~^, ■ ■ ■ , Xm^~^)'^)z- 

Now, TTi (resp. 712) can be chosen as the canonical projections from E"^ (resp. E™" — > 

jgm-i-j^ To be precise, if jc = (xi, X2, . . . , Xm)'^, then Tri{x) = xi, 7r2(a;) = (x2, • • • , Xm)^ ■ At 
this point, we are already able to compute the window set W numerically: Since the left PF~ 
eigenvector of the substitution matrix contains the the lengths of the prototiles (cf. remark 
after Thm. [TT^ above) . the lengths can be chosen such that V C Z[A]. That is, every x £V 
can be expressed as a; = ao + aiA + • • • + am-iX"^'^ , where Oi G Z. Then 

X* = (ao + 01^2 H + am-iX2^~^, ■■■ , ao + aiAm-i H h Om-iXm-i'^'^)'^ ■ 

More general, we may apply the star map to V, obtaining the (countable) set V* , which in 
turn yields the window set W = cl{V*). Actually, some of the pictures in the next section 
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were created in essentially in this way: Using the substitution to generate a large finite subset 
F C V, where the points are expressed as x above, we plot all points of F*. If this set F* 
appears to be bounded, even if the number of points is increased, then this is an indication 
that the considered set y is a model set. 

To prove rigorously that V is actually a model set, one proceeds by showing fi{W) > and 
IJ,{dW) = 0. If is a polyhedron, then this is an easy task. But in the generic case W is of 
fractal nature, and it may be even hard to prove that W is of positive measure. 

The last thing to prove is F = {■ki{x) \ x £ A,tt2{x) G W}, and not just a too small subset. 
More precisely: Let U := {7ri(a;) | x G A, 7r2(a;) G W}. Then it must hold dens(y) = dens(C/). 
In this context the following theorem turns out to be helpful |Scj . 

Theorem 1.7. Let V be a model set given by Qj. Then 

dens{V) = fi{W)/detA, 

where detA denotes the determinant of the generator matrix of A. Equivalently, detA is the 
volume of any measurable fundamental domain of A. 

1.3. Rauzy Fractals. As mentioned above, 'Rauzy fractal' is essentially just another term 
for the window set W. The term 'Rauzy fractal' occurs in the context of dynamical systems 
and number theory, while the terms 'window' or 'window set' arose from the theory of qua- 
sicrystals. Originally, 'Rauzy fractal' denoted one special fractal of this kind, namely the one 
arising from the substitution 

cj(l) = 12, cj(2) = 13, cj(3) = 1, 

which is analyzed in Therefore, the term 'generalized Rauzy fractal' is sometimes used 
to emphasize the fact that one deals with the more general case. (Of course, this is just a 
symbolic substitution, but it is trivial to formulate the geometric substitution according to 
Def. 11.11 cf. comment after (jlj)). 

According to |PYFj . the definition of a Rauzy fractal is essentially as follows: if a is a 
substitution for m prototiles in E^, where the factor is a PV number of algebraic degree m, 
then the Rauzy fractal is just the window set 1^ as in ©, where one does not care about 
f^{W) > or ^(dW) = or Theorem 11.71 Note that therefore it is much more difficult to 
show that a given substitution point set is a model set than it is to compute the corresponding 
Rauzy fractal. 

Before we proceed, let us mention a theorem which is important in this context and which is 
necessary in the next section. 

Theorem 1.8 ( |HUTj ). Let (X, d) be a complete metric space and {/i, . . . , /„} a finite set of 
contractive maps (i.e., 3c < 1 : Vj;,y G X : d{f{x),f{y)) < cd{x,y)). Then there is a unique 
compact set K G X such that K = |jr=i fi{^)- 

This theorem is a consequence of the Banach Fixed Point Theorem, and is easily generalized 
to partitions of K: Let {/i, . . . , /„} be as above and Ij C {1, . . . ,n} for 1 < j < ^, then 
there is a unique tuple {Ki, . . . , K() of compact sets Kj, such that Kj = Uie/j fii-^ji) 
1 ^ J ^ ^ QKLlj ). In the next section, this theorem will help to show the compactness of the 
window set W. 
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2. Two Rauzy Fractals and the Duality Principle 

We will describe the duality principle by choosing a certain one-dimensional substitution, 
constructing the corresponding Rauzy fractal and the associated dual tiling, and then we will 
use the duality principle to show that the dual tiling is a model set. 

2.1. Some Rauzy Fractals. It is easy to see, that for any nonnegative integer nxn-matrix A 
there is a substitution tiling in with substitution matrix A. Actually, from any such matrix 
one can immediately deduce a substitution, since there are only few geometric restrictions for 
tilings in one dimension. (The same problem for higher dimensions becomes quite difficult!) 
Here we consider the following matrices for n > 2: 

Mn := imij)i<ij<n, where rriij = 

Obviously all these matrices are primitive. A substitution with substitution matrix Mn for 
the prototiles 1, 2, . . . , n is given by 

1 ^ n 

2 n— 1 n 

n 1 2 • • • n— 1 n 

Note that this is just a symbolic substitution. But in one dimension the geometric realization 
is obtained from this scheme in a unique way: The prototiles are intervals of different lengths, 
and the lengths are given by the PF-eigenvector of M„. There is a nice trick to deduce the 
substitution factor directly from this substitutions: For fixed n G Z, n > 2, let ■= sin( ^^^j^ ). 
Because of the wonderful formula (see |NDj ) 

k-l 

(5) —Si = Si+l-fc+2i/ 

i/=0 

it follows = Sn, T^S2 = Sn-i + Sn, ■ ■ ■ , T^s„ = Si + S2 + - • • + Sn- Therefore the substitution 
factor (which is the PF-eigenvalue of M„) is j^, and the length of tile i is Sj. Now, we 
determine which of these substitutions define a Rauzy-fractal. The first condition to be 
fulfilled is that the substitution factor must be a PV-number^, cf. Theorem 11.61 This is the 
case for n = 2, 3, 4 and 7, and for no other n < 30. (This has been checked by a Maple 
program, and we don't expect any more PV-numbers as PF-eigenvalues of Mn for larger 
values of n.) The corresponding characteristic polynomials of M„ are: 

x^-x-1, rc^-2x^-x + l, x'^ -2x^ -3x'^ + X + 1, x"^ - 4x'^ - 62;^ + lOx^ + - - x + 1. 

The latter two are reducible over Z. The irreducible polynomials defining — (n = 2), ^ (n = 
3), M (n = 4) and ^ (n = 7) are 

(6) x^ -x-1, x^ -2x^ -x + l, x^ - 3x^ + 1, x^ - 4x^ - 4x^ + x + 1. 

Therefore, the algebraic degree of these numbers is 2,3,3,4, resp. Here we will focus mainly 
on the third case (n = 4). The first case (n = 2) leads to the well known Fibonacci sequences 



J 1 : i + j > n 
1 : i + j < n 



'Or a Salem number, but these don't play a role here, so we don't mention them explicitly. 
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(see Section 121 or |PYFj . Sections 2.6 and 5.4). The fourth case (n = 7) will lead to a three- 
dimensional Rauzy fractal in the end. The third case (n = 4) will show the methods and 
the typical problems which may occur, and therefore we discuss this case in the following in 
detail. This part may be regarded as a recipe how to compute other Rauzy fractals explicitly. 

2.2. The case n = 4. In this case we have got four prototiles. Since we consider one- 
dimensional tilings here, the tiles can be chosen to be intervals. We denote the four prototiles 
by L ('large'), M ('medium'), S ('small') and X ('extra small'). Their lengths are already 
known: They are given by a left PV-eigenvector of the substitution matrix, or by Q, so the 
values can be chosen as f^, f^, f^, = 1, resp. (or any positive multiples of these values). In 
order to construct the Rauzy fractal, we need this values as expressions in A = (cf- Section 
ll.2|l . W.l.o.g, we choose X := [0, 1]. The substitution now reads 

(7) X ^ L, S ^ ML, M SML, L XSML. 

Since AX = L it follows L = [0,A]. Let s,m,l denote the lengths of S,M,L, resp. Since 
^ = A, we conclude from ((Tj) the equations 

As = m + A 

Am = s + m + A = s(A + l) 
A^ = i + s + m + X = s(A + !) + ! = Am + 1 

and therefore m = ^^^^ . From © follows A^ = 3A^ — 1, hence A^^ = 3A — A^, and we obtain 

m = (A^ - 1) • (3A - A^) = A^ - 2A. 

Consequently, 

s = A2-A-m-l = A-l. 
To apply the standard construction of Rauzy fractals, the number of prototiles must equal 
the algebraic degree of the substitution factor. So we need to get rid of one letter. By ((Tj), a 
tile X in the tiling is always followed by a tile S. Since x + s = A, we can merge every pair 
XS into one tile L. It is easy to check that the new substitution 

(8) S ML, M SML, L LML 

gives rise to the same tilings as ^ does, up to merging pairs XS into L. We identify these 
tilings with a point set V, where V contains all (right) endpoints of the tiles, £ V, then 
V C Z[A], and every point x £ V can be expressed by x = a + 6A + cA^ {a,b,c E Z). The 
lattice A G now is given by 

(9) A = ((1, 1, 1)^, (A, A2, X^f, (A2, Xl Xlfh, 

where A2, A3 are the algebraic conjugates of A. From AA2A3 = — 1 and the fact that the values 
Sk occur frequently above, one easily deduces A2 = — t",A3 = This allows us to write 
down the star map now (cf. (dJ). Let £ V. If x £ V, then x = a + bX + cX^ for some 
a,b,c £ Z, and 

,10) ..^<.(;)o(^^).c(;;^:). 

By (jnj, the window set W is just the closure oiV*. A short computation (using Af = 3 A? — 1 
for i = 2,3) shows 
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In order to show the compactness of W, and possibly to draw further conclusions about W 
(which is the Rauzy fractal wanted), we consider a substitution set V such that a{V) = V. 
(As mentioned above, the considered substitution a can be easily formulated for the point set 
V as well as for the tiling T.) For example, let 



(11) 



V = ■■■ SMLLMLLMLSMLLMLMLSMLLMLLMLSMLL ■ ■ ■ 



where the right endpoint of the underlined L is located at zero. (Note, that the condition 
a{V) = V together with ((TT|) defines V uniquely.) Denote by Vl (W/j^s, resp.) the set of 
all points in V which are right endpoints of intervals L (M, S, resp.). Then obviously 



V = VlU Vm U Vs. 



Moreover, from 
(12) 



and cr{V) = V follows 



Vl 

Vm 

Vs 



XV 

XV- 

XVm 



XVl - A2 + a. 



A2 + A. 



This can be seen as follows: One tile of type M and one tile of type L occur in every 
substituted tile T, not depending on the type of T. If T = [a,b], then cr(T) contains tiles 
[A6 — A — m, Xb — A] (of type M) and [A6 — A, A6] (of type L). Therefore, from b G V follows 
Xb — X G Vm, yielding the second equation, and Xb G Vl, yielding the first part of the first 
equation. The other parts of (fT2|) are obtained analogously. 

Applying the star map to (|12() . we obtain the following equations for V*: 



V* 



VI 



u 



V* 



u 



V* 



(13) 



v, 



M 



VS 



QV* 
QV* - 

QV*M 



U QVl 

(A2,A3)^ 



(A2^A32)^ + (A2,A3)^^ 



(A2^A32)^ + (A2,A3)^ 



Obviously, Q defines a contracting linear map, thus all maps appearing in (|13|) are contracting. 
Therefore, Theorem 11.81 applies: There is one unique compact W = Wl U Wm U Wg which 
fulfills 1)131) . This assures the existence of this set, so we obtained a Rauzy fractal W . 





Figure 2. Left: The Rauzy fractal W arising from the substitution Q. Mid- 
dle: the same set, partitioned according to the distinction of tiles S, M and 
L. Right: The Rauzy fractal arising from the substitution Q for n = 3: 
S ^ L, M ^ ML, L SML. 
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To visualize W and its fractal nature, W can be approximated numerically. In fact, we just 
produce a finite subset of V and apply the star map to it. The resulting points are plotted 
into a diagram, which is shown in Fig. [21 (left). If we distinguish the points in the diagram 
according to their type {S,M or L), the result is shown in Fig. [2 (middle). The set W is the 
union of the three sets Wl,Wm,Ws, where Wi = cl(y*). From (|TS)) follows that Ws is an 
affine image of Wm, but not similar to Wm, and that Wm is an affine image of W. 

In the same way one can compute the Rauzy fractal for the case n = 3 in Q). The resulting 
pattern is shown in Fig. [21 right. Note, that we did not claim that one of the original 
substitution sets is a model set until now. In order to prove this, one needs to prove 
H{W) > 0, n{dW) = and dens{V) = dens({7ri(a;) | x G A,7r2(ic) G W}). 

2.3. The dual tiling. Recall that Q in (|13|) is a contraction. Therefore Q' := is an 
expansion. Applying Q' to l\l'S\i and replacing V* with its closure Wi {i G {S,M,L}) yields 
the following equation system: 

Q'Wl = WLUWMUWsUWL-{X2A3f + 

(14) Q'Wm = I^L - (1, 1)'^ U VFm - (1, 1)^ U l^s - (1, 1)"^ 
Q'Ws = Wm - iX2, XsV + {1, If 

This defines the substitution a' in with three prototiles Wl, Wm, Ws- 

(^'{{Wl}) = {WL,WM,Ws,WL-{X2,X3f + 

(15) (j'^Wm}) = {WL-{l,lf,WM-{l,lf,Ws-{l,lf} 

^'{{Ws}) = {WM-{X2,X3V+{l,m 

In contrary to Def. II. H the expansion by Q' is not a similarity, but an affinity. Anyway, these 
substitution rules will yield a tiling of E^. A part of it is shown in Fig. [21 (To be precise, 
we need to show that there will be no overlaps and no gaps, but this is not essential for 
the following, since we will consider the corresponding substitution point set, or a modified 
version where all tiles are polygons.) 




Figure 3. Left: A patch of the tihng arising from ^TEji . Tiles of type Ws 
are black, the other tiles come in two slightly different shadings, such that the 
boundary between tiles of the same type remains visible. Right: A polygonal- 
ized version. Tiles and colours correspond to the tiling on the left. 



Now, we will prove that the dual tiling is a model set. In analogy to the construction of V 
out of the original substitution tiling arising from we choose in every prototile a control 
point. A good choice is (cf. ((T?^): := G Wl, cm '■= Qcl — {X2,X3)^ = (A2,A3)"^ G 
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Wm, cs := Qcm — {^2^, ^3^)^ + (^2,^3)-^ G Ws- It is not too hard to check that the set 
{Wl, Wm -(A2^ X3^)^+2{X2, A3)^ + (l, 1)^} gives rise to a point set V such that a'{V') = V , 
and V is a Delone set. To be precise, let Vq := {cl, cm — (A2^, As^)-^ + 2(A2, As)-^ + (1, l)"^}, 
and a' be defined for points as well as for tiles in the canonical way. Then Vq C cj'(Vo'), and 
inductively {(t')''{V^) C (ct')''+H^o)- Therefore, 

(16) V'=[j{ar{K) 

k>l 

is well defined, and it holds (y'iV') = V' . 

Remark: As a first idea, one would try to start with the set {cl} instead of Vq. In this case, 
the property {cr')^{{cL}) C {cr')^~^^{{cL}) is fulfilled, but the resulting point set would be no 
Delone set, since it possesses arbitrary large holes. 

By the construction, it holds V C ((1, 1)^, (A2, As)^, (A2^, Aa^)^)^. In particular, this shows 
a curious property of V: If one coordinate of some z = [x, y)'^ G V' is known, say, x = 
a + 6A2 + cA2^ (a, b,c E Z), then the other coordinate y is known, too: y = a + bX^ + cAs^. 
Another unusual property is that each prototile occurs in just one orientation. This is not 
true for most known examples of nonperiodic substitution tilings in dimension larger than 
one. 

In a similar way as above (cf. (|12)) ) we obtain from (|15)) the following equations, where 
V' = VgUVlU VI and V- is the set of all points of type q in V {i £ {S, M, L}). 

Vi = Q'ViuQ'Vi,UQ'Vi-iX2,X3f + {l,lf 

(17) Vi, = QX-(A2,A3)^ug'V^-(A2,A3)^UQ'y^ 

n = Q'vi + {X2,X3f-2{X2^XsY^Q'vi, + {X2,X3V-2{X2^X3Y 

In order to show that V is a model set, we have to determine the lattice A' and the window 
set W for V' . By the construction of V follows A' = A. Since E V, we know that every 
cc E y is of the form a(l, 1)-^ + 6(A2, Aa)^ + c(A2^, As^)^. Thus our new star map is just the 
inverse of the original one in H1U|) . namely 

* : 7r2(A) ->■ E\ X* = 7ri(7r2"^(ic)) = a + bX + cA^. 

Applying -k to (|T7j) yields: 

v[* = x-^v[* u x-^vj^/ u x-^vi* - X + 1 

(18) VI/ = X-^V[* - XuX-^Vj^j* - XUX'^V^* 

= A-iy^* + A - 2A2 u A- V^/ + A - 2A2 

Since all occurring maps are contractions. Theorem 11.81 applies, and the unique compact 
solution is 

d{VL*) = [-A, 0] = L, cl{Vi/) = [A - A^, -A] = M, cl(yi') = [-2X^ + 1, -2X^ + X] = S. 

The multiplication of (|18j) by A yields a distinct geometric realization of our original substi- 
tution Q. So, the dual tiling of the dual tiling is the one we started with. In particular, the 
window set oi V is W = S U M U L = [-2X'^ + 1, -2A2 + A] U [A - A^, 0]. Therefore W is 
compact, 

(19) ^i{W') = A^ - 1 = 7.290859374 . . . > 

and fi{dW') = 0. All we are left with is to compute dens(y) and use Theorem 1 1.71 If the tiles 
Wi of the substitution a' in (|15|) would be polygons, this would be an easy task, using the 
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remark after Theorem EH Therefore we replace the fractal tiles Ws, Wm, Wl by appropriate 
polygons P5, Pa/, Pl. Let 

(A22-l,A32-l)^, b 

\2 — A2,A3 — As)-^, d 



a 
c 
e 
9 



(A2^ — A2, ) 
{2X2^ - A2 



1,2A32 



(A2,A3)J 



h 



{2X2' 



Xo 



I As' 



X, 
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1,2A32 



1 



(A2 + A2 — 1, A3 + A3 



11 



h-b 






Figure 4. Left: The construction of the polygonal prototiles. Middle: A 
sketch of the substitution for them. Right: The corresponding substitution 
matrix Af^i. Note that A^ji = {A^^)^ . 

Let conv(a;i, . . . , Xn) denote the convex hull of a^i, . . . , x^- Then (see Fig. |1} 



(20) 



Ps 
Pm 
Pl 



conv a 



g,c,e,e - g), 



b, a,h,h — b), 
conv(e — g,e,d,d — g) U conv(c 
conv(0, e, c, b) U conv(0, g, h, d, e) 

It is quite simple (but lengthy) to show that by replacing Wi with Pi in (|15|) . the substitution 
a' leads to a tiling, without any gaps or overlaps. (Essentially, one uses the definition of a 
substitution tiling in Definition 11.11 and the primitivity of the substitution matrix. To show 
that there are no gaps and no overlaps it suffices to consider all vertex configurations in the 
first n substitutions cr'{T), (a')'^{T), . . . , (cj')"(T) for appropriate n and some prototile T. So, 
in any certain case this is a finite problem. Here n = 4 and T = Pl will do.) It is also quite 
simple (but lengthy again) to compute the areas of the prototiles Ps, Pm, Pl- Unfortunately, 
we cannot use the fact that the left eigenvector of the corresponding substitution matrix 
reflects the ratio of the areas, since the substitution for the polygonal tiles does not preserve 
the shape any longer. E.g., cr'{Ws) is congruent to Q'Ws, but cr'{Ps) is not congruent to 
Q'Ps- Altogether, we obtain 



£1 

Si 



£2 _^ £1 

Si S2 



(211 



IJl{Ps) = I A2 1 A3 (A3 - A2) 
//(Pl) = A2" - A3' + 2IA2I + 2A3 - IA2IA3' - A3A2' = 1 + 



f) — A3 — A2 



£2 _^ £1 

Si S2' 



S4S1 

S2^ 



where the last equality requires an excessive use of ©• The vector {fi{Ps) , fJ,{PM) , h{Pl)) 
indeed is a left eigenvector of the substitution matrix A (cf. Fig. 0J) corresponding to the PF 
eigenvalue A. A right eigenvector corresponding to A is (,82, S3, 54)"^. Therefore, the normalized 
right eigenvector of ^ — which entries, by Theorem ll.31 are the relative frequencies pi of tile 
types Pi {i G {S,M,L}) in the tiling — is {ps,Pm,PlV '■= ^^+^+37(52, S3, S4). Since every 
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point of V corresponds to exactly one tile in the polygonal tiling, by Prop. 11.41 the relative 
density of V is 

dens(y') = {{^,{Ps),^l{PM),^^{PL)) ■ {ps,PM,PL f y^ 

(22) ( , , ^(S1S2^ . sj S2S3 S1S3 SiS4^\"^ 

^^"^^ = [S2 + S3 + S4) — 9 h — H \ h S4 H 5- 

V S4 S4 S4 S2 S2^ ) 

= 0.8100954858... 

In order to apply Theorem 1 1 . 71 we need to compute det A (cf. Q). A rather lengthy compu- 
tation (or the usage of a computer algebra system, here: Maple) shows det A = pTj^^^- 
Using elementary trigonometric formulas yields 

S1S2S4 = -(S3 + S4 - S2 - Slj 

1 / /37r\ \ 1 V3 

= -^S3 + 2cOS ^-J Si-SiJ=-S3 = — , 

and therefore det A = 9. Putting together this with (|19|) and (|22|) in view of Theorem 1 1 . 71 we 
obtain 

/i(W"')/dens(F') = 9 = det(A). 
The leftmost equality requires again a very lengthy computation, or the usage of Maple. This 
completes the proof, and we obtain the following theorem. 

Theorem 2.1. The set V in H^). arising from the substitution a' in is a model set. 

One may ask why we do not state that the original set V arising from the substitution (jH} is 
a model set. The reason is that there is one gap, namely, the value of ^^iW). For example, 
one needs to show that the tiling by the fractal tiles W5 , Wm , is not overlapping. This 
seems to be clear from the pictures, but we did not give a rigorous proof here. Indeed, this 
will require some more work, and we refer to further publications. Anyway, there are good 
reasons to assume that the areas of the fractal tiles Ws, Wm, Wl are equal to the areas of the 
polygonal tiles Ps, Pm, Pl, resp. Proving this will complete the proof that the original set V 
is a model set, too. 



3. Further Remarks 



The question whether a given structure is a model set is motivated for example by the question 
for the diffraction property of this structure. It is known that model sets have a pure point 
diffraction image, without any continuous part. For details, cf. BMRS . 

We did not emphasize it throughout the text, but all tilings considered here are nonperiodic. 
For the one-dimensional tilings this is a consequence of the fact that the relative frequencies 
of the tiles are irrational numbers. If a one-dimensional tiling is periodic, then the relative 
frequencies are rational numbers. Anyway, all occurring two-dimensional tilings are nonperi- 
odic, too. This can be shown with standard tools from the theory of nonperiodic substitution 
tilings (cf. statement 10.1.1 in |GSHj ). 

An important point in proving Theorem l2.1l is that we know the occurring values in terms of 
trigonometric expressions, thus it is possible to make exact computations of ^{W'), dens(y), 
det(A) etc. In general, the fractal nature of the occurring window sets may forbid this. But 
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if all occurring sets are polytopes the things become easier. In the following, some examples 
are listed where the duality principle provides short and simple proofs that these considered 
substitution tilings are model sets. 

Fibonacci sequences: The symbolic substitution S ^ L, L ^ LS gives rise to the 
widely examined Fibonacci sequences. (Not to be confused with the Fibonacci numbers 
1,1,2,3,5,8,13,. . ., sometimes also called Fibonacci sequence. If one considers the first substi- 
tutions of S, then the relation becomes clear: L, LS, LLS, LSLLS, LSLLSLSL, . . .). 

A geometric representation as a substitution tiling (Fibonacci tiling) is for example given by 
L := [0, 1], 5 = [1, r] (r the golden ratio) and 

a{{S}) = {L + T}, a{{L})={L,S}. 

By the standard construction (cf. Section fl. 2 j) . we obtain for some V which fulfills cr{V) = V 
(namely, V = (J^>g (T^'^({L, L — 1}), up to one point this will do) the equations 

VI = -T-^Vl U -T^'^Vs 

Vs = -T-^Vl + 1 
to be read as equations in E^. The unique compact solution is 

z\{Vl) = Wl = [-l,r-\ c\{V^) = Ws = [t-\t\. 
The dual tiling is given by the following substitution: 

<y'{{WL}) = {WL,Ws} 
<y'{{Ws}) = {WL-T} 

Note that the substitution factor is — r. Nevertheless, the resulting family of substitution 
tilings for a' is equal to these for a. That means the dual tilings of Fibonacci tilings are 
Fibonacci tilings again; or shortly: Fibonacci tilings are self-dual. Moreover, a short com- 
putation yields that the equation in Theorem 11.71 is fulfilled, and therefore the Fibonacci 
sequences are model sets. 



A 



A' 



A 




Figure 5. Left: The substitution for the Ammann chair. Right: A part of a 
corresponding substitution tiling. 

Golden triangles and Ammann chair: The substitution for the golden triangles (cf. Fig.^ 
defines a model set, if one replaces every triangle with the vertex at its right angle. This was 
essentially shown in |DvO) in detail^. It turns out that the golden triangle tilings are dual to 
the so-called Ammann chair tilings (or 'A2 tilings', cf. |GSHj ^. A sketch of the substitution 
of the latter is shown in Fig. [21 together with a part of the tiling. The substitution factor and 



The authors used all vertices in the tiling instead of one vertex of each triangle, but their result is easily 
transferred to the latter case. 
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the substitution matrix are equal to those of the golden triangle substitution. All occurring 
tiles, and all occurring window sets are polygons. Therefore, it is simple to compute the areas 
and the densities, thus proving — with the help of Theorem 11.71 — that both substitution 
tilings (resp. the corresponding point sets) are model sets. 

The duality of the Ammann chair tilings and the golden triangle tilings was already indicated 
in [(tELI briefly. But a proof of this needs some care. For example, one has to pay attention 
to the fact, that each prototile occurs in four different orientations, not just in one, as it is 
true for the tilings in Fig. IHI 

Variants of the tilings in Section |5J The tilings considered in the last section are a rich 
source of other substitutions, or other Rauzy fractals. If we restrict to the case n = 4, then the 
underlying lattice A will always be the same (namely, as in ©)• One possible modification 
is to permute the letters in the substitution. E.g., the substitution 



instead of ^ yields the Rauzy fractal shown in Fig. El (left). Another variant is to use the 
matrix A^i of Fig. \^ (the transpose of the original substitution matrix A^j) as a substitution 
matrix for a one-dimensional substitution with substitution factor A. One possibility is the 
substitution 



where S,M,L are intervals of length ^^,^2,54, resp. This substitution yields the Rauzy 
fractal shown in Fig. |^ (right). Note that in this case the set Wm is not connected, it 
consists of two parts which are similar to each other. The substitution tiling arising from 
this Rauzy fractal (resp. from the corresponding equation system, multiplied by Q' , cf. (|14j) ) 
would possess a prototile, namely, Wm, which is disconnected. Note, that we did not exclude 
disconnected (proto-)tiles in our definition of (substitution) tilings. 




Figure 6. Two more Rauzy fractals, the left one arising from the substitution 
H23|) . the right one arising from the substitution (|24|) 

Combining these two modifications leads to even more examples, all of them based on the 
lattice A. All the two-dimensional tilings are model sets, which is easily proven by using the 
results of Section [2 From the case n = 3 one can also construct different tilings by permuting 
the letters. In contrary, all variants of the case n = 2 will yield Fibonacci sequences. Anyway, 
in the cases n = 3, n = 4 it may be interesting to examine the different possibilities with 
respect to common properties of the corresponding Rauzy fractals. 



(23) 



S ML, M MLS, L LML 



(24) 



S ^ M, M ^ SML, L SMLL, 
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Substitution factor: rj = ^/t, 
where r = 



